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Abstract. Extensive photometric multisite campaigns of the 5 Scuti variable FG Vir are presented. For the years 
2003 and 2004, 926 hours of photometry at the millimag precision level were obtained. The combinations with 
earlier campaigns lead to excellent frequency resolution and high signal/noise. A multifrequency analysis yields 
79 frequencies. This represents a new record for this type of star. The modes discovered earlier were confirmed. 
Pulsation occurs over a wide frequency band from 5.7 to 44.3 c/d with amplitudes of 0.2 mmag or larger. Within 
this wide band the frequencies are not distributed at random, but tend to cluster in groups. A similar feature is seen 
in the power spectrum of the residuals after 79 frequencies are prewhitened. This indicates that many additional 
modes are excited. The interpretation is supported by a histogram of the photometric amplitudes, which shows 
an increase of modes with small amplitudes. The old question of the 'missing modes' may be answered now: the 
large number of detected frequencies as well as the large number of additional frequencies suggested by the power 
spectrum of the residuals confirms the theoretical prediction of a large number of excited modes. 
FG Vir shows a number of frequency combinations of the dominant mode at 12.7162 c/d (m = 0) with other 
modes of relatively high photometric amplitudes. The amplitudes of the frequency sums are higher than those of 
the differences. A second mode (20.2878 c/d) also shows combinations. This mode of azimuthal order m — -1 is 
coupled with two other modes of m — +1. 
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1. Introduction 

The 8 Scuti variables are stars of spectral type A and 
F in the main-sequence or immediate post-main-sequence 
stage of evolution. They generally pulsate with a large 
number of simultaneously excited radial and nonradial 
modes, which makes them well-suited for asteroseismolog- 
ical studies. The photometric amplitudes of the dominant 
modes in the typical S Scuti star are a few millimag. It is 
now possible for ground-based telescopes to detect a large 
number of simultaneously excited modes with submillimag 
amplitudes in stars other than the Sun (e.g., Breger et al. 
2002, Frandsen et al. 2001). Because photometric stud- 
ies measure the integrated light across the stellar surface, 
they can detect low-degree modes only. This is a simplifi- 
cation for the interpretation because of fewer possibilities 
in mode identification. 

A typical multisite photometric campaign allows the 
discovery of about five to ten frequencies of pulsation 
from about 200 to 300 hours of high-precision photometry 



(e. g., V351 Ori, Ripepi et al. 2003; V534 Tau, Li et al. 
2004) . These excellent observational studies are then com- 
pared to theoretical pulsation models, but the fit is hardly 
unique (e.g., 9 2 Tau, Breger et al. 2002b). The uniqueness 
problem can be lessened by studies with even lower noise 
in the power spectrum. This can be achieved by the very 
accurate measurements from space and by larger ground- 
based studies with more data, which concentrate on a sin- 
gle selected star. These more extensive ground-based stud- 
ies also lead to to higher frequency resolution. The latter 
is important, because 6 Scuti stars can show a large num- 
ber of very close frequency pairs (or groups), which can 
only be resolved through long-term studies lasting many 
months or years. The question of frequency resolution is 
an important aspect in planning asteroseismological space 
missions (e.g., see Handler 2004, Garrido & Poretti 2004). 

The Delta Scuti Network (DSN) is a network of tele- 
scopes situated on different continents. The collaboration 
reduces the effects of regular daytime observing gaps. The 
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network is engaged in a long-term program (1000+ hours 
of observation, 10+ years, photometry and spectroscopy) 
to determine the structure and nature of the multiple fre- 
quencies of selected S Scuti stars situated in different parts 
of the classical instability strip. The star FG Vir is the 
present main long-term target of the network. This 7500K 
star (Mantegazza, Poretti & Bossi 1994) is at the end of 
its main-sequence evolution. The projected rotational ve- 
locity is very small (21.3 ± 1.0 km s _1 , Mittermayer & 
Weiss 2003, see also Mantegazza & Poretti 2002). 

A number of photometric studies of the variability of 
FG Vir are available: a lower-accuracy study by Dawson 
from the years 1985 and 1986 (Dawson, Breger & Lopez 
de Coca 1995, data not used here), high-accuracy stud- 
ies from 1992 (Mantegazza, Poretti & Bossi 1994), as 
well as previous campaigns by the Delta Scuti Network 
in 1993, 1995 and 2002 (Breger et al. 1995, 1998, 2004). 
Furthermore, for the year 1996, additional uvby photom- 
etry is available (Viskum et al. 1998). 12 nights of data 
were of high accuracy and could be included. 

Because of the large scope of the long-term project on 
the pulsation of FG Vir, the photometric, spectroscopic 
and pulsation-model results cannot be presented in one 
paper. Here we present the extensive new photometric 
data from 2003 and 2004 as well as multifrcqucncy analy- 
ses to extract the multiple frequencies excited in FG Vir. 
The analyses concentrate on the available three years of 
extensive coverage (2002-2004) and also consider the pre- 
vious data (1992-1996). 

Separate studies, presently in progress, will (i) present 
mode identifications based mainly on high-dispersion line- 
profile analyses and the data presented in this paper, (ii) 
examine the nature of close frequencies, (iii) compute as- 
teroseismological models of stellar structure to fit the ob- 
served frequency spectrum. 

2. New photometric measurements 

During 2003 and 2004, photometric measurements of the 
star FG Vir were scheduled for ~ 350 nights at four ob- 
servatories. Of these, 218 nights were of high photometric 
quality at the millimag level with no instrumental prob- 
lems. These are listed in Table 1 together with the addi- 
tional details: 

1. The APT measurements were obtained with the T6 
0.75 m Vienna Automatic Photoelectric Telescope 
(APT), situated at Washington Camp in Arizona 
(Strassmeier et al. 1997, Breger & Hiesberger 1999). 
The telescope has been used before for several lengthy 
campaigns of the Delta Scuti Network, which con- 
firmed the long-term stability and millimag precision 
of the APT photometry. 

2. The OSN measurements were obtained with the 0.90 m 
telescope located at 2900m above sea level in the 
South-East of Spain at the Observatorio de Sierra 
Nevada in Granada, Spain. The telescope is equipped 
with a simultaneous four-channel photometer (uvby 



Stromgren photoelectric photometer). The observers 
for 2003 were: E. Rodriguez, P. Lopez de Coca, A. 
Rolland, and V. Costa. 

3. The SAAO measurements were made with the 
Modular Photometer attached to the 0.5 m and the 
UCT photometer attached to the 0.75 m telescopes of 
the South African Astronomical Observatory. The ob- 
servers were V. Antoci, E. Guggenberger, G. Handler 
and B. Ngwato. 

4. The 0.6-m reflector at Siding Spring Observatory, 
Australia, was used with a PMT detector. The ob- 
servers were P. Lenz and R. R. Shobbrook. 

The measurements were made with Stromgren v and 
y filters. Since telescopes and photometers at different ob- 
servatories have different zero-points, the measurements 
need to be adjusted. This was done by zeroing the average 
magnitude of FG Vir from each site and later readjusting 
the zero-points by using the final multifrequency solution. 
The shifts were in the submillimag range. We also checked 
for potential differences in the effective wavelength at dif- 
ferent observatories by computing and comparing the am- 
plitudes of the dominant mode. No problems were found. 

The measurements of FG Vir were alternated with 
those of two comparison stars. Details on the three-star 
technique can be found in Breger (1993). We used the 
same comparison stars as during the previous DSN cam- 
paigns of FG Vir, viz., CI = HD 106952 (F8V) and C2 = 
HD 105912 (F5V). No variability of these comparison stars 
was found. The two comparison stars also make it possible 
to check the precision of the different observing sites. The 
residuals from the assumed constancy were quite similar, 
i.e., for the (C1-C2) difference we find a standard devi- 
ation of ± 3 mmag for all observatories and passbands 
except for ± 2 mmag (2004 SAA075 v as well as y pass- 
bands) and ± 4 mmag (2004 APT75 v and 2003 OSN90 y 
measurements. The power spectrum of the C1-C2 differ- 
ences does not reveal any statistically significant peaks. 

The resulting light curves of FG Vir are shown in Figs. 
1 and 2, where the observations are also compared with 
the fit to be derived in the next section. 

3. Multiple frequency analysis 

The pulsation frequency analyses were per- 
formed with a package of computer programs 
with single-frequency and multiple-frequency 
techniques (PERIOD04 , Lenz & Breger 2005; 
http://www.astro.univie.ac.at/~dsn/dsn/Period04/), 
which utilize Fourier as well as multiple-least-squares al- 
gorithms. The latter technique fits up to several hundreds 
of simultaneous sinusoidal variations in the magnitude 
domain and does not rely on sequential prewhitening. 
The amplitudes and phases of all modes/frequencies 
are determined by minimizing the residuals between the 
measurements and the fit. The frequencies can also be 
improved at the same time. 
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Fig. 1. Multisite photoelectric three-star-photometry of FG Vir obtained during the 2003 and 2004 DSN campaigns. 
Ay and Aw are the observed magnitude differences (variable - comparison stars) normalized to zero in the narrowband 
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Fig. 2. Multisite photoelectric three-star-photometry of FG Vir obtained during the 2003 and 2004 DSN campaigns, 
continued 
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The file of the figure is too large to be uploaded to astro-ph. 



Please download the complete paper from 



http://www,a5tro.univie.ac.at/-dsn/ds5n/news2005/fgvir.pdf 



Fig. 3. Example of power spectra of the 1992-2004 data. Top: Spectral windows showing effects of the daily and 
annual aliasing. Bottom: New frequencies detected in the most difficult frequency region with the lowest amplitudes. 
The diagram shows that in the 40-45 c/d region pulsation modes are present and have been detected. The choice of 
which peaks are statistically significant depends somewhat on the details of the analysis. 



Our analysis consists of two parts: We first examine 
the extensive 2002 - 2004 data and then add the available 
1992 - 1996 data. 

3.1. Frequencies detected in the 2002 - 2004 data 

The following approach was used in an iterative way: 

(i) The data were divided into two data sets to sepa- 
rate the y and v filters, each covering the total time pe- 
riod from 2002 - 2004. This is necessary because the am- 
plitudes and phasing of the pulsation are strongly wave- 
length dependent. In principle, the different amplitudes 
could be compensated for by multiplying the v data by an 
experimentally determined factor of 0.70 and increasing 
the weights of the scaled v data accordingly. (Anticipating 
the results presented later in Table 2, we note that this 
ratio is confirmed by the average amplitude ratio of the 
eight modes with highest amplitudes.) However, the small 
phase shifts of a few degrees cannot be neglected for the 
larger-amplitude modes. Consequently, the data were an- 
alyzed together for exploratory analysis, but not for the 
final analyses. 

(ii) We started with the single-frequency solution for 
the two data sets using the program PERIOD04. For the 
Fourier analyses the two data sets were combined to de- 
crease the noise, while for the actual fits to the data, sep- 
arate solutions were made. 



(iii) A Fourier analysis was performed to look for 
additional frequencies/modes from the combined residu- 
als of the previous solutions. Additional frequencies were 
then identified and their signal/noise ratio calculated. 
Following Breger et al. (1993), a significance criterion of 
amplitude signal/noise = 4.0 (which corresponds to power 
signal/noise of ~12.6) was adopted for non-combination 
frequencies. The most clearly detected additional frequen- 
cies were included in a new multifrequency solution. In 
order for a new frequency to be accepted as real, the 
signal/noise criterion also had to be fulfilled in the mul- 
tifrequency solution. This avoids false detections due to 
spill-over effects because the Fourier technique is a single- 
frequency technique. Furthermore, since there exist reg- 
ular annual gaps, trial annual alias values (separated by 
0.0027 c/d) were also examined. We note that the choice 
of an incorrect annual alias value usually has little or no 
effect on the subsequent search for other frequencies. The 
choice of an incorrect daily alias (separated by 1 c/d) 
would be more serious and we carefully examined different 
frequency values. 

(iv) The previous step was repeated adding further fre- 
quencies until no significant frequencies were found. Note 
that only the Fourier analyses assume prewhitcning; the 
multiple-frequency solutions do not. 

In this paper we omit the presentation of very lengthy 
diagrams showing the sequential detection of new frequen- 
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cies, except for the example shown in the next subsection. 
A detailed presentation of our approach and its results can 
be found in our analysis of the 2002 data (Breger et al. 
2004). 

FG Vir contains one dominant frequency: 12.7162 c/d 
with a photometric amplitude five times higher than that 
of the next strongest mode. To avoid potential problems 
caused by even small amplitude variability, for this fre- 
quency we calculated amplitudes on an annual basis. The 
results of the multifrcquency analysis are shown in Table 2. 
The numbering scheme of the frequencies corresponds to 
the order of detection, i.e., the amplitude signal/noise ra- 
tio, and therefore differs from that used in previous papers 
on FG Vir. 

3.2. Further frequencies detected in the 1992 - 2004 
data 

An extensive photometric data set covering 13 years is es- 
sentially unique in the study of S Scuti stars, promising 
new limits in frequency resolution and noise reduction in 
Fourier space. The noise reduction is especially visible at 
high frequencies, where the effects of systematic observa- 
tional errors are small. The analysis of the combined data 
had to work around two problems: the earlier data is not 
as extensive as the 2002 - 2004 data and there exists a 
large time gap between 1996 and 2002. 

The time gap did lead to occasional uniqueness prob- 
lems for the frequencies with amplitudes in the 0.2 mmag 
range: next to the annual aliasing of 0.0027 c/d we find 
peaks spaced 0.00026 c/d, corresponding to a ~ 10 year 
spacing (see Fig. 3, top right). Fortunately, the excellent 
coverage from 2002 - 2004 minimized these ambiguities. 

The relatively short coverage of the data from 1992 and 
1996 excluded the computation of 79-frcqucncy solutions 
for individual years (to avoid overinterpretation) . We have 
consequently combined all the y measurements from 1992 
to 1996 as well as the available 1995 and 1996 v data. 
Together with the y and v data sets from 2002 - 2004, we 
had four data sets. 

Fig. 3 illustrates some examples of the resulting power 
spectra. Due to the large amount of APT data from 2002 
- 2004, which therefore dominates, the 1 c/d aliases are 
not zero (top left). Nevertheless, due to the excellent fre- 
quency resolution, these aliases are very narrow so that 
the aliasing problem is not severe. The figure also shows 
the power spectrum of the measurements in the 40 - 45 
c/d region, which was the most difficult region for us to 
analyze due to the small amplitudes of all the detected 
frequencies. 

The new detections are included in Table 2. 

A comparison with the frequencies published in ear- 
lier papers shows that all the previously detected frequen- 
cies were confirmed. This also includes those previously 
detected modes not found to be statistically significant 
in the 2002 data alone. In a few cases, different annual 
aliases were selected. However, the main result is the in- 



crease in the number of detected frequencies to 79, which 
more than doubles the previous results. 

Fig. 4 shows the distribution of frequencies in fre- 
quency space. We note the wide range of excited frequen- 
cies, which is unusual for 5 Scuti stars, as well as the clus- 
tering of the excited frequencies. This clustering persists 
even after the suspected combination frequencies and 2f 
harmonics are removed. 

A new feature is the detection of frequencies with val- 
ues between 40 and 45 c/d. They all have small amplitudes 
of Ay ~0.2 mmag. The lower noise of the new data now 
made their detection possible. 

3.3. Color effects 

The light curves of pulsating stars are not identical at dif- 
ferent wavelengths. In fact, amplitude ratios and phase 
shifts provide a tool for the identification of nonradial 
modes (e.g., see Garrido et al. 1990, Moya et al. 2004). 
For 6 Scuti stars, the amplitude ratios between different 
colors are primarily dependent on the surface tempera- 
ture. For the individual pulsation modes, the phase dif- 
ferences and deviations from the average amplitude ratio 
are small. This means that observational errors need to 
be small and any systematic errors between the different 
colors should be avoided. 

For most nights there exist both v and y passband 
data, so that amplitude ratios as well as phase differences 
can be derived. However, our 79-frequency solution is not 
perfect. In order not to introduce systematic errors in the 
phase differences and amplitude ratios, for the calculation 
of amplitude ratios and phase differences, we have omitted 
those nights for which two-color data are not available. 
Consequently, no data from 1992 and 1993 were used and 
all 1995 (single-color) CCD measurements were omitted. 

Table 3 lists the derived phase differences and am- 
plitude ratios for the modes with relatively high ampli- 
tudes. The uncertainties listed were derived from error- 
propagation calculations based on the standard formulae 
given by Breger et al. (1999). The results can now be used 
together with spectroscopic line-profile analyses to iden- 
tify the pulsation modes. 

4. Combination frequencies 

We have written a simple program to test which of the 79 
frequencies found can be expressed as the sum or differ- 
ences of other frequencies. Due to the excellent frequency 
resolution of the 2002 - 2004 data, we could be very restric- 
tive in the identification of these combinations. A generous 
limit of ± 0.001 c/d was adopted. The probability of incor- 
rect identifications is correspondingly small. A number of 
combinations was found and these are marked in Table 2. 
They generally agreed to ± 0.0002 c/d. 

How many accidental agreements do we expect? We 
have calculated this number through a large number of 
numerical simulations, assuming a reasonable agreement 
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Table 2. Frequencies of FG Vir from 2002 to 2004 



Frequency Detection' 1 - 1 Amplitude Frequency Detection Amplitude 

cd _1 Name Type Amplitude v filter y filter cd _1 Name Type Amplitude v filter y filter 

S/N ratio millimag S/N ratio millimag 











± 0.04 


±0.05 










± 0.04 


±0.05 


5.749f 


f47 




5.7 


0.48 


0.28 


23.3974 


fl2 


(2) 


22 


1.73 


1.25 


7.9942 


144 




5.9 


0.42 


0.34 


23.4034 


u 


(2) 


71 


5.65 


4.02 


8.3353 


178 


fe-fi 


3.8 


0.33 


0.17 


23.4258 


Us 




4.9 


0.37 


0.28 


9.1991 


h 




53 


3.82 


2.78 


23.4389 


f24 




9.0 


0.68 


0.52 


9.6563 


f 5 




71 


5.20 


3.61 


23.8074 


£48 




5.5 


0.45 


0.32 


10.1687 


f 26 




8.6 


0.64 


0.42 


24.0040 


f 5 8 




4.6 


0.35 


0.31 


10.6872 


179 


f 4 -fi 


3.5 


0.26 


0.15 


24.1940 


fio 




29 


2.26 


1.58 


11.1034 


f 2 (> 




11 


0.67 


0.65 


24.2280 


f 3 




74 


5.79 


4.20 


11.2098 


f38 




6.4 


0.35 


0.43 


24.3485 


fl8 




12 


0.99 


0.63 


11.5117 


fro 


fs-fi 


4.0 


0.30 


0.18 


24.8703 


f36 


fl+f 2 


6.4 


0.48 


0.38 


11.6114 


fco 




5.3 


0.35 


0.27 


25.1788 


f30 




7.3 


0.59 


0.39 


11.7016 


f33 




6.9 


0.42 


0.44 


25.3793 


f 45 




5.7 


0.37 


0.31 


11.8755 


f63 




4.4 


0.27 


0.25 


25.4324 


fl5 


2fi 


16 


1.23 


0.89 


11.9421 






7.6 


0.55 


0.38 


25.6387 


f68 




4.0 


0.30 


0.26 


12.1541 


h 


(2) 


85 


6.09 


4.21 


26.5266 


in 




4.7 


0.30 


0.20 


12.1619 


fl4 


(2) 


16 


1.13 


0.83 


26.8929 


fei 




4.5 


0.34 


0.27 


12.2158 


f 6 2 




5.1 


0.38 


0.24 


26.9094 


in 




4.2 


0.34 


0.19 


12.7162 


fl 




442 


31.74 


21.92 


28.1359 


fl9 




11 


0.87 


0.55 


12.7944 


fl7 




13 


0.88 


0.66 


29.4869 


fsi 


f 7 +fn 


5.2 


0.39 


0.28 


13.2365 


f27 




8.3 


0.45 


0.56 


30.9146 


f60 




4.5 


0.31 


0.27 


14.7354 


f49 




5.3 


0.29 


0.38 


31.1955 


f 5 7 




4.6 


0.35 


0.27 


16.0711 


fl3 




20 


1.56 


1.07 


31.9307 


f 3 4 




6.6 


0.49 


0.40 


16.0909 


fsi 




7.3 


0.55 


0.39 


32.1895 


£32 




7.0 


0.56 


0.38 


19.1642 


f 26 




8.6 


0.55 


0.59 


33.0437 


f 7 4 


(4) 


4.3 


0.29 


0.19 


19.2278 


f 9 




30 


2.51 


1.69 


33.7677 


f 43 


fi+fs 


5.9 


0.44 


0.31 


19.3259 


f41 




6.3 


0.53 


0.34 


34.1151 


f22 




9.8 


0.75 


0.49 


1 A C A o n 

19.6439 


ffi5 




4.3 


0.40 


0.22 


34.1192 


r 

172 




4-7 


0.21 


0.25 


19.8679 


fs 


(3) 


55 


4.44 


3.19 


34.1864 


f 5 4 




4.9 


0.37 


0.25 


19.8680 




(3) 


30 


2.40 


1.78 


34.3946 


f 5 5 




4.6 


0.31 


0.28 


20.2878 


fll 




26 


2.13 


1.45 


34.5737 


f 23 




9.3 


0.69 


0.45 


20.2925 


fr>6 




4.6 


0.31 


0.39 


35.8858 


f 76 




4.1 


0.22 


0.21 


20.5112 


f 35 




6.6 


0.41 


0.52 


36.1196 


£40 


fl+f 4 


6.3 


0.40 


0.31 


20.8348 


f39 




6.3 


0.53 


0.38 


36.9442 


£37 


fl+f 3 


6.4 


0.43 


0.27 


21.0515 


f 6 




55 


4.43 


3.08 


39.2165 


f69 




4.0 


0.27 


0.16 


21.2323 


fl6 




14 


1.06 


0.80 


39.5156 


f59 


f 9 +fll 


4.5 


0.24 


0.25 


21.4004 


f46 




5.7 


0.52 


0.33 


42.1030 


f64 


2f 6 


4.3 


0.22 


0.28 


21.5507 


f28 




7.9 


0.61 


0.42 


42.1094 


f62 




4.5 


0.24 


0.25 


22.3725 


f42 




6.2 


0.52 


0.35 


43.0134 


f73 




4.4 


0.21 


0.15 


23.0253 


f66 




4.1 


0.30 


0.23 


43.9651 


£S7 




4.0 


0.26 


0.17 


23.3943 


f 2 l 




11 


0.85 


0.60 


44.2591 


f77 




4.0 


0.20 


0.18 


(1) The noise 


for the amplitude 


signa 


1/noise ratios were calculated 


over a 4 cd _1 


range. 


Limits for 


a significant detection 


are 



4.0 for independent frequencies and 3.5 for combination modes with known values. Numbers in italics indicate 1992 - 2004 
data (see text). 

' 2 ' The close frequencies, 12.1541 and 12.1619 as well as 23.3974 and 23.4034 c/d, are all separate modes. In short data sets 
this could lead to an erroneous identification as single modes with variable amplitude. 

For the possible frequency pair near 19.868 c/d the existence of two separate modes cannot be proved at this point. A 
single frequency with a slowly variable amplitude (beat period ^21.5 years) is also possible. 

^ The 2002-2004 data clearly show a mode at 33.044 c/d, though with considerably reduced amplitudes from 1992-1995 
data. Breger et al. (1998) listed the value of the frequency as 33.056 c/d, which was the highest peak from a broad selection of 
peaks separated by annual aliases 0.0027 c/d apart. We note that in the new data, a value separated by 1 annual alias, viz., 
33.0461 c/d, is also possible. 
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Fig. 4. Distribution of the frequencies of the detected modes. The diagram suggests that the excited pulsation modes 
are not equally distributed in frequency. 

Table 3. Phase differences and amplitude ratios 



Frequency Phase differences in degrees 

C/d (j) V — (j>y 

2002-2004 1995-2004 1 ' 



Amplitude ratios 
v/y 

2002-2004 1995-2004 



A 


12.716 


-1.7 ± 0.1 


-1.5 


1.45 ± 0.00 


1.45 


h 


12.154 


+3.1 ± 0.7 


+2.8 


1.44 ± 0.02 


1.44 


h 


24.227 


-3.0 ± 0.7 


-3.0 


1.38 ± 0.02 


1.40 


A 


23.403 


-3.5 ± 0.7 


-3.3 


1.40 ± 0.02 


1.42 


A 


9.656 


-5.1 ± 0.8 


-4.5 


1.44 ± 0.02 


1.43 


h 


21.051 


-3.7 ± 1.0 


-4.2 


1.44 ± 0.02 


1.45 


A 


9.199 


-6.7 ± 1.1 


-6.6 


1.38 ± 0.03 


1.40 


h 


19.868 2) 


-4.3 ± 1.9 


-2.9 


1.45 ± 0.05 


1.44 


h 


19.228 


-4.6 ± 1.7 


-3.9 


1.48 ± 0.05 


1.46 


ho 


24.194 


-1.2 ± 1.9 


-0.6 


1.43 ± 0.07 


1.38 



1 ' Error estimates omitted: usually lower than for 2002-2004, but some instability is possible due to large time gap 
2 ' Using single frequency with annual amplitude variations 



of the observed frequency to within 0.0002 c/d of the pre- 
dicted frequency. We obtain an average of 0.93 acciden- 
tal matchings of peaks with combination frequencies. We 
conclude that most or all detected combination frequen- 
cies are not accidental. The argument is strengthened by 
the fact that the combinations detected by us all contain 
one of two specific modes, which reduces the chance of 
accidental agreements to essentially zero. 



We also note that the lowest frequency detected, f 4 7 at 
5.749 c/d, can be expressed as a triple combination of /i, 
f3 and f57. This may be accidental. 

4.1. Combinations of the dominant mode at 12.7163 

c/d 

Due to the presence of a dominant mode at 12.7163 c/d 
(fi), it is not surprising that some combination frequen- 
cies, fi±fi, exist and are detected (see Table 2). In order 
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Fig. 5. Distribution of the amplitudes of the frequencies with significant detections. To increase the accuracy, we have 
computed amplitudes from 0.5*(y amplitude + 0.70 v amplitude) to simulate the amplitude in the y passband. Note 
the large number of detected modes with amplitudes near the detection limit of 0.2 mmag. This suggests that even 
moderate increases in the amount of data lead to considerably higher number of detections. 
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Fig. 6. Power of the residual noise of the 1992 - 2004 data after subtracting 80 (79) frequencies. (The additional 
frequency is the close doublet at 19.868 c/d adopted for the long time span.) The average amplitude in the amplitude 
spectrum was calculated for 2 c/d regions and squared to give power. The strong increase towards lower frequencies 
is caused by observational noise and standard in terrestrial photometry. Note the excess power in 10 - 50 c/d region, 
shown above the dotted curve. This shows that many additional pulsation modes exist in the same frequency region 
in which the already detected frequencies occur. 



to examine this further, we have performed additional cal- 
culations with the 2002 - 2004 data. We have repeated our 
multifrequency solutions described earlier while omitting 
all frequency combinations of fi . The residuals of the y and 
v data were combined to form the sum (y + 0.70v) to ac- 



count for the different amplitudes at the two passbands. A 
new multifrequency solution containing the possible com- 
binations of the dominant mode with f 2 through f 8 was 
made. The results are shown in Table 4, in which small 
differences compared to Table 2 are caused by the different 
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Table 4. Combination frequencies involving the dominant 
mode 



fl± 


Amplitude 


in y (2002-2004) 




Sum of frequencies 


Difference of frequencies 




mmag 


mmag 


h 


0.35 


(0.06) 


h 


0.28 


0.19 


U 


0.31 


0.17 


f 5 


0.33 


(0.16) 


fe 


0.39 


0.21 


h 


(0.10) 


(0.07) 


U 


(0.07) 


(0.13) 



The amplitudes in brackets are too low for fulfilling the 
adopted criterion of statistical significance. 



procedures of our analyses. The amplitudes of the sums 
(fi+f») are higher than those of the differences (fi-fj). We 
believe the result to be real and intrinsic to the star. The 
differences are generally found at low frequencies, where 
the observational noise is higher. Increased noise should 
lead to higher amplitudes, which are not found. 

4.2. Combinations of the 20.2878 c/d mode 

Apart from the mode combinations involving the dom- 
inant mode, fi, two other two- mode combinations are 
found, both involving fn at 20.2878 c/d. While the mode 
identifications are still in progress, this mode can to a high 
probability be identified as a £=1, m=-l mode (Zima et 
al. 2003). At first sight, this appears surprising, since the 
photometric amplitude is only 1.5 mmag. However, for the 
known inclination and mode identification, we calculate a 
geometric cancellation factor of ~3.6, so that the real am- 
plitude of this mode is 5 mmag or larger. 

The fact that the m value is not equal to zero has some 
interesting consequences for the observations of combina- 
tion frequencies. The reason is that there are two frames 
of reference: the stellar frame corotating with the star, and 
that of the observer. For nonradial modes of m values ^ 
(i.e., waves traveling around the star), the frequencies 
between the two frames of reference differ by mfi, where Q 
is the rotation frequency of the star (see Cox 1984 for an 
excellent discussion). The frequency combinations occur 
in the corotating (stellar) frame, and not the observer's 
frame of reference. Consequently, many possible combina- 
tions involving non-axisymmetric modes should not be ob- 
served as simple sums or differences of observed frequency 
values. 

It follows that for the non-axisymmetric (m=-l) mode 
at 20.2878 c/d, our simple method to search for combina- 
tion frequencies from the observed frequency values may 
only detect combinations, fj + fj, with m = +1 modes. 

This strict requirement is met by the two identified 
combinations of 20.2878 c/d! Both coupled modes, fg = 
19.2278 c/d as well as fn = 20.2878 c/d have been iden- 



tified as m = +1 modes. Such a combination of m values 
of opposite sign can be detected because the combination 
is invariant to the transformation between the two frames 
of reference: — mfl + mfl = 0. 

5. The problem of missing frequencies solved? 

5 Scuti star models predict pulsational instability in many 
radial and nonradial modes. The observed number of low- 
degree modes is much lower than the predicted number. 
The problem of mode selection is most severe for post- 
main-sequence S Scuti stars, which comprise about 40 
percent of the observed S Scuti stars. The theoretical fre- 
quency spectrum of unstable modes is very dense. Most 
modes are of mixed character: they behave like p-modes 
in the envelope and like g-modes in the interior. For ex- 
ample, for a model of the relatively evolved star 4 CVn, 
the models predict 554 unstable modes of I — to 2, i.e., 

6 for i = 0, 168 for t = 1, and 380 for i = 2 (see Breger 
& Pamyatnykh 2002). However, only 18 (and additional 
16 combination frequencies) were observed (Breger et al. 
1999). The problem also exists for other 5 Scuti stars. A 
complication occurs since the models so far cannot predict 
the amplitudes of the excited modes. 

Two explanations offer themselves: the missing modes 
exist, but have amplitudes too small to have been ob- 
served, or there exists a mode selection mechanism, which 
needs to be examined in more detail. Promising scenarios 
involve the selective excitation of modes trapped in the 
envelope (Dziembowski & Krolikowska 1990) or random 
mode selection. 

Let us turn to the star FG Vir. Unpublished mod- 
els computed by A. A. Pamyatnykh predict 80 unstable 
modes with £ = 0, 1 and 2 in the 8 - 40 c/d range. This 
number is smaller than that mentioned previously for the 
more evolved star 4 CVn, but until now this large number 
was not observed either. The present study addresses the 
question by lowering the observational amplitude thresh- 
old to 0.2 mmag. We have detected 79 frequencies, of 
which 12 could be identified as harmonics or combina- 
tion frequencies. This leaves 67 independent frequencies. 
There also exists considerable evidence that many more 
modes are excited: 

(i) Consider the amplitude distribution of the detected 
modes shown in Fig. 5. There is a rapid increase in the 
number of modes as one goes towards low amplitudes. 
The present limit near 0.2 mmag is purely observational. 
Consequently, the number of excited modes must be much 
larger. 

(ii) Consider the power spectrum of the residuals after 
subtraction of the multifrequency solution (Fig. 6). We see 
excess power in the 10 - 50 c/d range. This is exactly the 
region in which the previously detected modes were found. 
This indicates that many additional modes similar to the 
ones detected previously are excited at small amplitude. 

We can exclude the possibility that the large number 
of observed frequencies is erroneous because of imperfect 
prewhitening due to amplitude variability. We have ex- 
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amined this possibility in great detail, with literally thou- 
sands of different multifrcquency solutions allowing for 
amplitude variability. In no case was it possible to sig- 
nificantly reduce the structure in the power spectrum of 
the residuals. In fact, the 'best' multifrcquency solution 
adopted treated the two colors as well as the 1992 - 1996 
and 2002 - 2004 data separately and allowed annual am- 
plitude variability of the dominant mode at 12.7162 c/d. 
Amplitude variability, therefore, cannot explain the excess 
power. 

Consequently, the problem that the number of detected 
modes is much smaller than the number of predicted low-t 
modes no longer exists, at least for FG Vir. Of course, we 
cannot conclude that each theoretically predicted mode 
is really excited and has been detected. This would re- 
quire much more extensive mode identifications than are 
available at this stage. 

In the previous discussion, we have concentrated on 
the low-£ modes which are easily observed photometri- 
cally. One also has to consider that at low amplitudes, 
variability from modes of higher £ values might also be 
seen photometrically. The geometric cancellation effects 
caused by the integration over the whole surface only be- 
come important for I > 3. This is shown by Daszyhska- 
Daszkiewicz et al. (2002), who calculated the amplitude 
reduction factors caused by temperature variations across 
the disk. From this paper we can roughly estimate a can- 
cellation factor of ^50 in the y passband for £ = 3 modes, 
implying that only the largest-amplitude modes could be 
photometrically detected. Such modes would have ampli- 
tudes similar to, or larger, than that of the observed dom- 
inant {£ = 1) mode at 12.7162 c/d and might therefore 
be expected to be few. Regrettably, the situation is some- 
what more complicated. The results presented in Fig. 2 of 
the Daszyhska-Daszkiewicz et al. paper are actually based 
on models with higher surface temperatures and could fit 
the (3 Cep variables. A. Pamyatnykh has kindly calculated 
specific models fitting the star FG Vir. Here the geomet- 
ric cancellation factor becomes smaller by a factor of two 
or three. Consequently, some of the low-amplitude modes 
observed by us could also be £ = 3 modes. 

We conclude that the large number of detected fre- 
quencies as well as the large number of additional fre- 
quencies suggested by the power spectrum of the residu- 
als confirms the theoretical prediction of a large number of 
excited modes. A mode-by-mode check for each predicted 
mode is not possible at this stage. 
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